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AN ITERATIVE METHOD DERIVED FROM EXISTENCE AND UNIQUENESS 

THEOREMS FOR SYSTEMS OF SECOND-ORDER, NONLINEAR, 
TWO-POINT-BOUNDARY-VALUE DIFFERENTIAL EQUATIONS 

By 

R. L. Manicke 
Manned Spacecraft Center 

s 

SUMMARY 

This paper develops some applications of the contrac- 
tion mapping principle* It shows sufficient conditions to 
guarantee existence of a unique solution of a two-point- 
boundary-value system of differential equations. Moreover, 
it exposes an iterative method that must converge from any 
initial guess to the unique solution, provided certain 
criteria are satisfied. An example is given to clarify the 
procedure . 

INTRODUCTION 

A boundary-value problem in Ordinary differential 
equations differs from an initial-value problem in that it 
must satisfy conditions at two or more points. An nth 
order differential equation can have n conditions at one 
point or one condition at n points. In fluid mechanics 
and trajectory mechanics , systems of second-order equations 
with conditions specified at two points are often encoun- 
tered. This paper will consider only systems of this form. 
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SYSTEM WITH IMPOSED CONDITIONS 


Consider 


Cl) 



yj(t) ♦ f x (t, y 1 (t), y 2 (t), 

y 2 (t) * f 2 (t, y i (t), y 2 (t), 

• % 

• • 

• • 

yj; it) ♦ f n (t, y^t), y 2 (t), 

* 1 ^ 


•••» y n (t), y{(t). yj(t), 
•••» y n (t), yJCt). yj(t) f 

•••» y n (t), yJCt), yj(t), 
• , y A (b) - 


y;(t» 


y;(t)) 


y^c*)) ■ 


o 


o 




0 


or in vector notation 

7 " + 7 (t, 7 » 7 ') * 0 with T ■ (y x » y 2 , *•*» y„) T 

s 

y"(a) = _7 , 7(b) = "b" 

When f -satisfies a Lipeohitz condition t, it means 
that each component satisfies the following: 


I f ± (t , y x . y 2 . •••» y^ •••> y n . yj. y{> •••* y„) 

- f i (t, y x , y 2 , •••, x it •••, y n , y[, •••» x’, •••, y^)| 


s K t |y ± - x.| + L.|y! - xj | with K., L i Z 
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and (t, y x , y 2> y if y n , y{, yj, y^) , 

(t, y x , y 2 , •••, •••» y„f y[* •••» •••, y^) elements 

of the domain of f i - [a,b] * R 2n , 1 S i S n . Define 

the minimum K , that will preserve the inequality as the 
Lipschitz constants of f^ . Notice that this condition 
does not presuppose the existence of a solution yl(t) of 
system (I).., 

GREEN’S FUNCTION FOR THE SYSTEM 

Now, in order to construct an iterative procedure, it 
is sufficient to use a Green’s function to transform each 
equation of a homogeneous boundary-value system into an 
equivalent integral equation system. (For the definition 
and properties of Green’s functions, see No. 2 in Biblio- 
graphy.) Such a Green’s function is: 

( & - - a < s s t s b 

G(t,s) = < 

( 0> :-|)iL- - a) - a £ t $ s S b 

(This Green’s function is for the operator -y”.) Since 
zero-boundary conditions are rarely encountered, it is 
usually necessary to transform the system into one with 
zero-boundary conditions and satisfy the original Lipschitz 
condition . Such a transformation is 

T(t) « alT - bX +_ (7 - ~B)t 
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Clearly, r i (a) - and T i (b) « . Therefore, if 7(t) 

is a solution of problem (I), then 7(t) ■ "y"( t ) * 7(t) is 
a solution of 


(II) 


7" + f (t , 7, 7' ) ■ o 
7(a) * 0 « 7(b) 


where 


T(t, T,T') ■ 7(t, 7 + 7(t) , 7' + 7' (t)) 

* 

Since 


F,(t, 


V 


T * 


* 


F, (t, 


U » 






^(t, •••• ♦ T^t), •••, I n ♦ T n (t), •••, I’ ♦ ♦ T n (t) 

- f t (t. •••• q t ♦ TjCt), •••, X n ♦ T n (t), •••, q' * - i (tl , t; * T n (t)| 


S K l J(i i ♦ t 4 (t)) - (q t ♦ t 1 (t))| ♦ t l |(*i * T'CD) - (q- ♦ T t (t))| 
* K i| l i ’ <»l| * L i| I i ‘ q i| 


this is exactly the same Lipschitz condition that each 
£ i , 1 £ i & n , must satisfy. 
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Theorem: 


Problem (II) is equivalent to 


/•b ^ 

T(t) • J G(t, s)F (s , z(s), z'(s)) ds 


( 1 ) 


(Two problems are equivalent if and only if they have 
exactly the same solutions.) 


Proof: 


7(t) 


f G(t, s)7(s, T(s) , T 1 (s)) ds 
•'a 

k f (s - a)7(s, T(s) , T* (s)) ds 
•'a 

+ g J* (b - s)7(s, T(s), T’(s)) ds 


Clearly, T(a) * 0 * 7(b) . 

By the product rule for differentiation, it may be 
found that 

7'(t) - b j 1 a J (s - a)7(s, T(s) , T' (s) ) ds 


(t - a)T(t, z(t) , z'(t)) 
>b 




s)F(s, z(s) , z' (s)) ds 


+ S— (b - t)7(t, 7(t), T'(t)) 


IT^T 
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ft 


7’(t) * (t ‘ *)7(t, 7(t), 7'(t)) 

+ ( - 2 1 ♦ a + b)7(t, T(t) , 7* (t)) 

♦ *t~: y (t - a)7*(t, 7(t) , 7» (t) ) 

+ ' 7(t , 7(t) , 7' (t) ) 

+ ( 2 1 - a - b)7(t , T(t) , 7* ( t ) ) 

♦ - ( a ^ ( r * ”*"(1) , 7* (t) ) 

■ - ]>1 7(t,7(t),7(t)) 

• -7(t.. 7(t), 7'(t)) 

or 

7"(t) * -F(t, 7(t) , 7' (t)) 

The fundamental theorem of contraction mappings will 
be applied to this class of systems. (A function 

T : S-iH£2*.s where S is a normed linear space is said 
to be a contraction mapping if and only if there is a real 
number ot,0 < a < 1 , such that for all x, y e S 
1 1 Tx - Tyll s o||x - y| | .) 
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CONTRACTION MAPPING THEOREM 


For every contraction mapping defined on a complete 
rtormed linear space S , there exists a unique fixed point 
x e S such that Tx * x . Moreover, for any x Q e S , 
the iterates 


T n x. 


T n-1 x 


T n “ 2 x 2 * ••• 


0 " ”1 

converge to the fixed point x , and 

.n 


Tx 


n- 1 


n 


X - X 
n 


T^S " x i * x o 


(See Bibliography No. 3, pp. 213-216.) 


THEOREMS FOR UNIQUE EXISTENCE 


First consider the system "y" + f(t,y) * 0 , ")T(a) * A, 
and !T(b) * T . Let S be the space of continuous functions 
on [a,b]. This linear space is complete with norm 


ax ) max | v (t) | I 
iSn \[a f b] x f 


max 
ISiSn 


In the following theorems it will be necessary to obtain 

b r b 


upper bounds for J G(t,s) ds and J [ f ' S ^ 

a a 


ds . 


Lemma 1: 


/. 


b G(t,s) ds s ( b 


- / a b | 2S % ii l * * 


b - a 
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Proof; 


A 




a) ds 



s) ds 


lb - t m - al l* (t »?(b - t) 2 

2 (b - a) 2 (b - a) 


% 


Cb - tUt - a) < (b - a) 2 
2 * S ' 


and 



3 G ( t , s ) 
5t 



a) 


ds + 




(b - s) ds 


(t - a) 2 + (b * t) 2 ^ b - a 

"2 (B ■-" ay * ~T~ 


Theorem 1 ; 

If the system "p* + 7(t,")T) * 0, y"(a) * A, and 7(b) * ”? 

satisfies the Lipschitz condition t and b - a < (8/K) 1//2 

where K « max K. , then a unique solution exists* 

IjJign x 

Proof : 

First, the system must be transformed into a homogeneous 
system; it then becomes *?" + 7(t,T) * 0 , Ti(a) * 0 « !T(b) . 
With the right side of equation (1) as the operator in the 
following equation, an attempt is made to find sufficient 
conditions when it is contracting. Notice that the operator 
T defined by 


8 


it ** 


t(f) * f G(t,s)F(s, T(s) , 7* (s)) ds is a linear mapping 

J a 


from S into S . 


r b 

Tu i (t) - Tv^Ct) - J G(t,s)[F i (s, u x (s), u 2 (s), •••, u^s), 


u n ( s )) ‘ F i( s » v i (s), 

v n C s )] ds 


Now, the Lipschitz condition is used to find 

/•b 

iTu^t) - Tv i (t)| S J G(t,s)K i |u i (s) - v^s)! 


_ __ /*b 

|Tu.(t) - Tv i (t) | £ Kjju - v|| J G(t,s) ds 


Since by Lemma 1 


f h 

J G(t,s) 


ds £ 


b - a 


therefore 


|Tu - T7| J < K| |u - T| | 


b - a 


, K * max k. 

l^iSn 1 


The mapping is contracting when 


K Y^? - a < 1 
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. •• J*- 



or if 


b - a < (8/K) 1/2 

a unique solution exists, namely: ylft) * T(t) + T(t) . 

Now, consider the system "y* M + f(t, y, y*) * 0, 

)T(a) * A, and )T(b) = B . This time, let S be the space 
of continuously differentiable functions. The space is 

complete with norm | |v| | * max 
Theorem 2 : 

Mp m "W^ m ___ 

If the system y" + f(t, y, y f ) * 0, y(a) « A, and 
Tib) * TT satisfies the Lipechitz condition t, and 


» A m 

with 

K = max K. 
lgign 1 

K(b - a) 2 + L(b - a) 
5 Z 

< i , 


and 

L ■ max L. 

l£i£n 1 


then a unique solution exists. 

Proof : 

The same transformation as Theorem 1 shall be used; 
again using T(t) , the system is transformed into a new one 
with homogeneous boundary conditions: 

7' + T(t, 7, 7 f ) = o 
7(a) = 0 = 7(b) 


max (k |v (t) I + L |v' (t) 
[a,b] V 1 1 1 x 
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Now form, 

m b 

|Tu i (t) - Tv i (t)| s J G(t,*)|f 1 (», u 1 (»), •••, u i (s), •••, « n C*). u[(s). «•*. u ^Cs)) 

- i L (s* v I (s), •••, v 1 (s), •••, v(n), v^(s), •••, v^(5))| ds 
$ /* b G(t > 5)CK 1 |u i (s) - v^s)! ♦ Lju’(s) - v’(s)|) ds 

s | par - Til / b G(t,») ds i I ITT - T| I 

•'a 


and 


n Tu i (t) • A Tv i(*) 


/ b |»^sl|| fi(s , ..., u . (s)) 


- ^(s, •••, V^(s)) 


ds 


Hence by Lemma 1 


which together imply 


r b i 

3G(t,s) i 

A 

at j 


S 14-2. | |T - T 


| | Tu - Tv| | S ^ a ^ 2 + L ^ - - g - 2) . 


|TT - v 


where K * max K. , L * max L. 

lgign 1 lgi^n 1 
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***** *41 ^ • ■ * m < 


Therefore, if 


K(b all + Mb - jj 


< 1 


the mapping ' T is contracting. Hence, a unique solution 
exists, namely: T( t ) * 7(t) + T(t) . 


ITERATIVE PROCEDURE 

If the system is of the form T 1 * + F(t,*z) * 0 , 

T(a) * 0 « 7(b) , the second part of the contraction mapping 

( 0 ? 

theorem states that for any starting value z(t) e S , 
the iterates of the ith component converge to z ± (t ) , that 
is 


(n) 


rb 

/ Gft 

•ra 


I (n-l) 

,s)F i ls, z 1 (s), 


(n-l) 


•••» z^s), 


(n-l) \ 

> 2 n (s >) 


ds 


CO) /( 0) (0) (0) \ T 

converges to z i (t) for any 7(t) =lz x (t), z 2 (t), •*•, z n ( t )l 
or 

(n) _/ (n^l)\ 

? + F 1 1 , z ) * 0 

(n) (n) 

T(a) = 0 = 7(b) 


converges to a solution. 
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EXAMPLE 

Consider the system 

(1) y’£(t) + y x (t) + y 2 (t) - o 

(2) y' 2 (t) + y x (t) + sin y 2 (t) = 0 

y x (0) = l y x (i) = l 

y 2 (0) = o y 2 (i) = o 


i 

i 


so that 


Ifift, y t , y 2 , •••, y t ' •••, y„, y x , yj) 

- £i(t, y x » y 2 , •••, x A . •••, y n . y ± . •••. y;)l 

S K.|y. - x. | 
is for (1) 

+ y 2 ■ x i - y 2 i * (1) iyi -'*J 

and for (2) 

I y i + sin y 2 - y 1 - sin x 2 | g (1) | y 2 - x 2 | 



T 2 (t) = 0 , SO ?(t) = (1,0) T 

The system is transformed into the homogeneous system 
z'^Ct) + z x (t) + z 2 (t) +1 = 0 

z"(t) + z r (t) + 1 + sin z 2 (t) = 0 

since ?(t, T, T*) = T(t, T + T(t) , z' +T(t)) 
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( 0 ) 

Let T * (1, tt/2) . Now, using the Green’s function 

for this system, it is found that 


( 1 ) 

z x (t) 


J f G(t,s)(2 + tt/2) ds 
0 

(1 - t) f s (2 + tt/2) ds + (- 1) f (s - 1) (2 + tt/2) ds 

•'0 


t 2 (l - t) (4 + tt) . t(s - l) 2 




Ml! 

M) 


i-p [t 2 - t 3 + t 3 - 2t 2 ♦ t] 


4 + TT 

3 


t(l - t) 


Similarly , 


CD 

z 2 (D 


_ 3t (1 - t) 


so 


(D 

T 


- t), 
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By the contraction mapping theorem, 




1/8 
- 1/S 


( 1 , */ 2 ) 




+ 

T 


IT 


t(l - t) 


3t(l - t) 

■MMM' 


r 



max {max 

[ 0 , 1 ] 




where i « 1,2 


2 .385 


Since this must converge uniformly, the next iterate obtained 

( 1 ) ( 0 ) 

from using T as a new T must be even "closer, 1 ’ (which 
is made precise by the definition of the norm). 


CONCLUSION 

This investigation has shown how an iterative method 
to solve two-point-boundary-value systems of differential 
equations is obtained from the principle of contraction 
mappings. Moreover, it states some conditions that will 
guarantee this solution to be unique, and that the method 
will uniformly converge to this unique solution for any 
starting value. Possible logical extensions of these 
results could be found by assuming a variable but bounded 
terminal value b . In cases where the length b-a is too 
large to guarantee unique existence, find conditions that 
can determine specific intervals on the real axis for which 
unique solutions exist. 
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